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,_^^ ' We perform a forecast of the CMSSM for the LHC based in an improved Bayesian analysis 

»^ ■ taking into account the present theoretical and experimental wisdom about the model'— '. In 

this way we obtain a map of the preferred regions of the CMSSM parameter space and show 
that fine-tuning penalization arises from the Bayesian analysis itself when the experimental 
value of Mz is considered. The results are remarkable stable when using different priors. 
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The start of the LHC has motivated a lot of effort to try to anticipate which kind of physics 



Q^' beyond the Standard Model is more likely to be there. Since the present experimental data are 



not powerful enough to select a small region of the parameter space of SUSY models, Bayesian 
0"| Statistics becomes a very powerful tool to try to make an inference of the probability of certain 

r^i ' regions of parameters of these models, where the choice of judicious prior probability for the 

parameters becomes more relevant. 
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10 ' 1 Bayesian Statistics 

(N' 

iVJ I The probability density of a particular point {p^} in the parameter space given a certain set of 

• ' data is the posterior probability density function (pdf), given by 

8 : ,(rf|dau) ^ ^i^^^i^. (1) 

where p(data|p^) is the likelihood, the probability density of measuring the given data for the 
K> ', chosen point in the parameter space. p{p^) is the prior, the "theoretical" probability density 

^ I that we assign a priory to the point in the parameter space, p(data) is the evidence. If one 

is interested in comparing regions of the parameter space of a given model, p(data) is just a 
normalization constant. 

2 Bayesian approach and Naturalness 

The parameters of the MSSM should not be far from the electroweak scale in order to avoid 
unnatural fine-tunings to obtain the correct scale of the electroweak breaking. From the mini- 
mization of the tree-level form of the scalar potential. 
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M'z = 2 ^\ "' , ' - 2^^l^ , (2) 



we can see that if fi and the soft masses rriH^ 2 ^^^ ^o^ close to the electroweak scale, a big 
cancellation is necessary in order to obtain the right value of Mz- A conventional measure of 
this cancellation are the fine-tuning parameters'-^—, 
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Since naturalness arguments are actually statistical, one may expect that a penalization of fine- 
tunings should arise from the Bayesian analysis itself. Let us see how this comes about. Let us 
consider Mz on a similar foot to the rest of experimental data, 



p(data|s, m, M, A, B, fi) = Cmz -Crest 



(4) 



where s represents the SM parameters, £rest is the likelihood of the all physical observables 
except Mz, and Cmz is the likelihood of Mz- Let us now use the sharpness of the likelihood 
of Mz to approximate Cmz — ^{Mz — M^^^) and marginalise the pdf in the ^u— parameter, 
performing a change of variable /i — > Mz '■ 
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where /xq is the one which reproduce the experimental value of Mz as a function of the rest of 
the parameters Q. Now, comparing this to the definition of fine-tuning parameters one gets. 



pis,m,M,A,B\ data) = 2— C^cst pis,m,M,A,B, jjo) 

Mz Cf, 



(6) 



where the presence of the fine-tuning parameter does indeed penalize regions of the parameter 
space with large fine-tunings. As we will see, this is enough to make the high-energy region of 
the parameter space statistically insignificant. 

We have performed a Bayesian analysis of the MSSM with the following improvements'^ 
as we show above, the fine-tuning penalization arises from the Bayesian analysis itself; we have 
made a rigorous treatment of the nuisance variables where Yukawa couplings are fundamental 
parameters in contrast with previous analysis; and, finally, we have used an efficient set of 
variables to scan the MSSM, {m, M, A, B, mu, yt} ^ {m, M, A, tan /?, Mz,mt}. The last change 
of variables introduces a jacobian factor. 



pigi,mt,m,M, A,tan(3\data) = C^cst J\^^=^lo Pi9i,yt,m, M,A,B, ^i = no) 
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In addition, we have developed sensible priors which assume that soft-breaking term share a com- 
mon origin. This analysis have been implemented using MultiNest"^ algorithm as implemented 
in the SuperB ay eS^ code which incorporate SoftSusy , SusyBSG , Superlso— , MicrOMEGAs 
'^ codes. 

As can be seen in Fig 1, besides making the high-energy parameter space quite irrelevant, 
the EW breaking has a another remarkable effect, the probability distribution (pdfs) based on 
logarithmic or flat prior are quite similar after the incorporation of the EW scale. 

3 Experimental Constraints 

In this section we incorporate all the relevant experimental information to the likelihood piece 
of the probability distribution. We start by considering the most reliable and robust pieces of 



^We have ignored here the nomahzation factor. See eq. [T] 
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Figure 1: ID marginalized posterior probability distribution of the M for logarithmic (left panels) and flat (right 
panel) priors in the fi > case, for a scan including the EW breaking information (A/^^^). The red lines represents 

the marginalized prior. 

experimental information: EW and B(D)-physics observables and lower bounds on the masses of 
super symmetric particles and the Higgs mass. The left panel of Fig. 2 shows the pdf of {M — m} 
plane once this experimental information is incorporated. Clearly the bulk of the probability is 
now pushed into higher energy. This effect is basically due to the Higgs mass bound. It is well 
known that the tree- level Higgs mass is bounded from above by Mz, so radiative corrections are 
needed. Concerning the other observables, everythin g wo rks fine as long as SUSY is not at too 
low scale. We also show the discovery reach of LHC^^ for 1 fb^^ and 100 fb^^. These lines 
correspond to A = and tan/3 = 45, but they provide a good indication for the LHC discovery 
potential. 

Next, we have added information about anomalous magnetic momentum of the muon, a^. 
Taking e~^e~ — > hadrons data, there is a 3.3a discrepancy between the experimental and SM 
theoretical prediction, which has been often claimed as a signal of new physics. If one accept 
this, the supersymmetric masses should be brought to quite small values in order to produce 
a large enough contribution J'^^^^a^ to reconcile theory and experiment, as we can see in the 
center panel of Fig. 2. If we instead use r— decay data there is no big discrepancy, SUSY 
contribution does not need to be large and the pdfs are essentially unchanged by its inclusion. 

Supersymmetry offers good candidates for Cold Dark Matter (CDM), the most popular and 
natural one is the lightest neutralino. The right panel of Fig. 2 shows the pdfs after the as- 
sumption of CDM is made of neutralinos and without the information from a^. 

We have repeated the same analysis for ^u < 0. We have evaluated the ratio of the evidence 
of both cases, the Bayes factor, in order to compare the relative probability of the /i > and 
/x < branches. The /i < is slightly favoured, due to its capability to reproduce the central 
value of & — )• S7, but the effect is not really significant. On the other hand a^ favours /z > 
branch, this effect is stronger when i^DM is included. This is because ^dm constraints favours 
the low-energy region of the parameter space, and this is strongly preferred by a^. 

In conclusion, LHC offers an exciting horizon for SUSY discovery, but there is still a possibil- 
ity that this escapes detection, especially if the Higgs mass is not close to its present experimental 
bound. 
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Figure 2: 2D marginalized posterior probability distribution for logarithmic (upper panel) and flat (lower panel) 

priors in the n > case including: EW + B(D)-physics observables (left panel); + a^ (center panel); + CDM 

(right panel). The inner and outer contours enclose respective 68% and 95% joint regions. 
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